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Abstract 

A vector-like colorless fermion doublet and a singlet added to the Standard Model 

allow a consistent interpretation of dark matter in terms of the lightest neutral 

- . particle, as they may help in obtaining successful gauge coupling unification. We 

^ ' analyze in detail the mass range of the lightest neutral particle below the W mass, 

(•~^ . i.e. in a range of the parameters where the physics of the Standard Model Higgs 

CsJ ' boson may be substantially affected either directly or indirectly. 

> . 
Q . 

O ! 1 Introduction 

^ : 

PsJ , The Standard Model (SM) of the electroweak interactions is more than 30 years old 

' and it has been able to reproduce with great precision the many experimental results 

■ obtained until now. In particular at LEP the theory was tested at the per mille level 
without finding any discrepancy with the theoretical predictions. However, in spite 
of this extraordinary success, we are convinced of new physics beyond the SM, since 
there are problems where the SM does not provide an adequate solution. One of these, 
supported by observations, is the lack of a dark matter candidate. 

\ The most direct and impressive evidence of the existence of dark matter are surely 

^ ' the flat rotation curves of spiral galaxies. Other evidences for dark matter were found 

. at different scales, from galactic scales (several kiloparsecs) and clusters of galaxies 

^ \ (Megaparsecs) to global scales of hundreds of Megaparsecs pp. The total matter density 

■ can be inferred from the measurements of the power spectrum of the fluctuations in 
ly^ . the Cosmic Microwave Background (CMB). The recent measurements of the Wilkinson 

\ Microwave Anisotropy Probe (WMAP) [2j have shown that the total matter abundance 

^ ' in the universe is VLmh^ = 0-1277l']QgQ79; these measurements have provided also the 

^ ■ baryon abundance which is Q.bh? = 0.02229 ± 0.00073. We conclude that all the matter 

• ^ ! in the universe cannot be baryonic and that the dark matter abundance is Qdm^'^ = 

X ■ 0.1054l°-°n?q- The fact that a significant part of dark matter must be non baryonic was 

■ known before WMAP measurements. Indeed an estimate of was already available 
[3], and the value of was inferred from primordial nucleosynthesis [5; the difference 
seemed notable also at that time. An additional evidence for the non baryonic nature 
of the dark matter is given by structure formations: in a universe with only baryons the 
primordial density perturbations have not had enough time to grow and generate the 
galaxies observed today in the sky. These observations, however, do not tell us anything 
about the particle nature of dark matter. Then the question is about the nature, the 
origin and the composition of this important component of our universe, since dark 
matter does not find an explanation in the framework of the Standard Model of particle 
physics. 

Particle physics provides us with a large number of dark matter candidates, which 
appear naturally in various frameworks for reasons completely independent from the 
dark matter problem, and certainly not invented for the sole purpose of explaining the 
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presence of dark matter in our universe. Among these candidates an important distinc- 
tion is between particles created thermally or not thermally in the early universe. For 
thermal relics another important distinction is about how they decoupled from the pri- 
mordial soup, in particular if they were relativistic (hot dark matter) or nonrelativistic 
(cold dark matter). Arguments from large structures make us believe that a large, and 
presumably dominant, fraction of dark matter is made of cold relics. A well motivated 
class of cold dark matter particles are the so called Weakly Interacting Massive Par- 
ticles (WIMPs), which have mass between 10 GeV and a few TeV and interact only 
through weak and gravitational interactions, because the limits on charged relics are 
very stringent [5j. 

Another missing opportunity for the Standard Model is that gauge couplings do not 
quite unify at high energy; a possible solution is to add weakly interacting particles to 
change the running, in order to make unification work better. 

In this work we discuss a model that has both a cold dark matter candidate and can 
improve considerably over the Standard Model in the direction of successful gauge cou- 
pling unification. We introduce new matter with respect to the Standard Model alone, 
and we restrict ourselves to the case in which the added particles are fermions. Adding 
just a vector doublet allows remarkable improvements for unification; this model, fur- 
thermore, is highly constrained since it contains only one new parameter, the Dirac 
mass for the degenerate doublets, whose neutral components are the dark matter candi- 
dates. Such model, however, is ruled out by direct detection experiments: the vector-like 
vertex with the Z boson for the neutral particles remains unsuppressed, giving a spin- 
independent cross section that is 2-3 orders of magnitude above current limits This 
drawback can be solved by including a fermion singlet, with Yukawa couplings with the 
doublets and the Higgs boson. Doing this we generate a mixing between doublets and 
singlet, so that the neutral particles become Majorana fermions which have suppressed 
vector-like couplings with the Z boson. We assume a parity symmetry that acts only on 
the new fields. This imposes that they do not couple to ordinary matter. It also implies 
that the lightest particle is stable and, if neutral, it constitutes a good dark matter can- 
didate. This model has been introduced in reference [7J where a detailed dark matter 
analysis for high values of the relic particle mass can be found. In reference [7j it is also 
shown how the gauge coupling unification at high energy can be achieved and a rate for 
the proton decay is predicted that could be tested in the future. 

In this work we focus on the region of parameter space where the mass of the lightest 
neutral particle (LNP) is smaller than the W boson mass niw ~ 80 GeV. The analysis 
for higher mass was already done, as said above, but the main reason for doing so is that 
well above the WW production threshold, in order to account for the entire dark matter 
abundance observed, the mass M of charged components of the doublets is quite high. 
An important fact is that, for relatively low values of the LNP mass, the effects on Higgs 
boson physics are significant, both direct and indirect. On the one hand there are new 
available decay channels for the Higgs boson, and the decays into neutral particles may 
dominate the total width. On the other hand the new particles contribute to electroweak 
observables, so that they may change the indirect upper limit on the Higgs mass and 
improve the naturalness of the Higgs potential [8]. There are thus reasons to give special 
attention to this region of parameter space. 

The structure is the following: in section [2] we present the model with its spectrum, 
in section [3] we compute the relic abundance of the dark matter candidate, in section H] 
we discuss direct detection, in section [5] the effects on Higgs boson physics. Finally in 
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section [6] we consider a possible CP violating phase giving rise to an electron electric 
dipole moment. Conclusions are given in section [71 



2 The model 

The model consists of the following extension of the Standard Model 

C = CsM + AC (1) 

where we add to the Standard Model lagrangian the following renormalizable lagrangian 
(other than the kinetic terms for the various new fields) 

AC = XFHS + X^FcH^S + MFF^ + fiS^ + h.c. (2) 

The doublets Fc and F have respectively hypercharge ±1/2, is a singlet and H is the 
Standard Model Higgs doublet. We introduce the symmetry 

F, F„ S ^ -F, -F„ -S (3) 

with all other fields invariant. This imposes that the new fields do not couple to ordinary 
matter. We suppose the parameters (A, A'^, M, fj,) to be real (in sectionElwe will consider 
the effects of introducing a phase). The physical fields are chosen as follows 

The components of Fc and F are left-handed Weyl fields. The Goldstones and x can 
be put to zero by choosing the unitary gauge. 

In the charged sector there is a simple Dirac term of mass M, hence we define the 
Dirac spinor ip = F^ + {F~Y 

In the neutral sector we define the fields Ni as 

iVi = i= (fO - fO) = (fO + fO) N, = S (5) 

so that the mass matrix takes the form 

Mn ={ -M -^/2av (6) 

\ -\/2/3u -yfiav -2fi I 

where the Yukawa couplings have been replaced by the parameters 

« = ^ /? = ^ (7) 

We have to find now eigenvalues and eigenvectors of this matrix; in the general case this 
task cannot be accomplished analytically and therefore we will diagonalize the mass 
matrix numerically. Let mj be the eigenvalues and let V be the matrix that performs 
the diagonalization. We define Xi as the eigenvector corresponding to mj, i.e. 

Ni = VijXj V^MnV = diag (mi, m2, mg) (8) 

We identify the lightest neutral particle (LNP) with the index Z, then xi is the field of 
the LNfIi]. 

^from now on the index I for xi indicates lightest, and it must not be confused with left 
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3 Dark matter analysis 



In this section we compute the thermal rehc abundance of the LNP using the standard 
formahsm [9j. Before proceeding we should justify why we can use it, because there are 
situations in which this method fails [lOj . We have checked that in the parameter region 
of interest to us the masses of the other two neutral and of the charged particles are 
far higher than the LNP mass itself, so we can neglect coannihilations. The standard 
method is also not valid when the relic particle lies near a mass threshold since the LNP 
particles are Boltzmann distributed. Given our LNP mass range the only threshold 
present is that for WW production. For rrii > 75 GeV the WW process suppresses the 
LNP relic abundance to an unacceptable level, whereas for mi < 75 GeV it can be safely 
neglected. 

The evolution of the LNP number density ni is governed by the Boltzmann equation 



— + 3Hni 
at 



n, 



(9) 



where H is the Hubble parameter, is the LNP equilibrium number density, v^ei 
is the relative velocity and < avrei > is the thermal average of the annihilation cross 
section. The relevant temperatures are of order mi/25, so the Boltzmann equilibrium 
distribution is well justified. The Boltzmann equation can be solved approximately. 
First we introduce the variable x = mi/T. Secondly we parameterize the temperature 
dependence of the annihilation cross section as 



< aVrel >=(TqX~ 



(10) 



where n = corresponds to s-wave annihilation, n = 1 to p-wave annihilation, etc. At 
early times ni is accurately approximated by n^'^ , but as the temperature drops below 
the mass m;, n^'' drops exponentially until a point denominated "freeze out" is reached 
where the reaction rate is not fast enough to maintain equilibrium. From this point 
on, the n^'^ term in equation Q can be neglected and the remaining equation is easily 
integrated. Thus the solution of ([9]) is given by solving in two regimes and matching 
those solutions at the freeze out. The value of the freeze out point Xf \s obtained by 
imposing the equality between the interaction rate F = niavrei and the expansion rate 
H, and it is given by the numerical solution of the following equation 



Xf + ln + 



1 



Inx 



f 



In 



0.038 (n + 1) [g/gl^"^] mpiniiao 



where the Planck mass is mpi = 1.22 x 10^^ GeV and g^: is the number of effectively 
relativistic degrees of freedom at the time of freeze out. The present mass density of the 
relic particles is expressed as 



(n + l) ^ 



Xf 0.034 pb 



1/2 

9* 



(12) 



In our model we are dealing with a cold relic, therefore in the early universe, just 
before the decoupling, thermal equilibrium is maintained via LNP annihilations into 
fermions. There are two possible processes: Z boson exchange and Higgs boson ex- 
change, both p-wave. Cross sections for these processes and their thermal averages are 
given in appendix [Al In the following discussion we fix the values of Yukawa couplings 
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and analyze the model as function of {fi, M) for each case. The limit of small Yukawa 
couplings A, A'^ is not interesting, since in this case the LNP coincides approximately 
with the singlet, and the only way to produce all the dark matter observed is with 
the LNP mass near the Z pole or the Higgs pole. If A*^ = A the model possesses a 
SU{2)l X SU{2)ji symmetry broken to SU{2)v by the Higgs vacuum expectation value, 
and the coupling with the Z boson is suppressed. Also in this case the only way to 
produce all the dark matter is near the Higgs pole. 

We consider for the complete analysis two limiting cases: almost equal Yukawa 
couplings (symmetric case, or more properly nearly symmetric) and when one of them 
is vanishingly small (asymmetric case). To be consistent with negative searches from 
LEP we assume rrii > 4:5GeV and M > 100 GeV. The cases which we discuss are (the 
reason for doing so is explained in section [5]) 

• Symmetric case: a = 1.0 and /3 = 0.1 

• Asymmetric case I: a = 0.5 and /? = 0.5 

• Asymmetric case II: a = 0.65 and /? = 0.65 

Before proceeding we must say something about the Higgs boson mass, since the 
annihilation cross section for Higgs exchange depends on it and we have to choose its 
value carefully. We will see in section [5] that in the symmetric case the corrections to 
the electroweak parameter T are negligible, then the indirect upper limit on the Higgs 
mass valid in the Standard Model {nih ^ 166 GeV^ at 95% CL [11|) remains unchanged. 
On the contrary, in the asymmetric cases T is strongly affected by the new particles, so 
the upper limit is raised. We choose the reference values 

• Symmetric case: nih = 120 GeV 

• Asymmetric case I and II: nih = 300 GeV 

We plot our results in the {fi,M) plane in figure [TJ In the symmetric case only the 
relative sign of /i and M is physical, and our convention is M > 0. In the asymmetric 
case both signs are unphysical, then we choose /i > and M > 0. We identify the 
parameter space region for which 45 GeV < mi < 75 GeV and inside it we shade the 
area for which 0.089 < Qih'^ < 0.122 (corresponding to the 95%CL region from WMAP 
[2]). In all the cases the dark matter abundance can be accounted for by our LNP. 

Another check must be done: neutral particles could have been produced at LEP2 

e+e~ XI Xni (13) 

where the index nl stands for "next to lightest". Given the assumed symmetry ^ the 
only allowed decay for Xni is 

Xnl^Xlfl (14) 

where / indicates a generic fermion and / the corresponding antifermion. Since no such 
event was seen this may constrain the model. We have checked that it has not been 
kinematically allowed at LEP2, since the "next to lightest" particle mass is always above 
200 GeV in the parameter space region of interest. 
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Figure 1: LNP relic abundance. Contours for mi = 45, 75GeV are denoted by the solid 
lines, shaded regions correspond to 0.089 < Qih^ < 0.122 (WMAP 95%CL region) . 

4 Direct detection 



Dark matter particles of the Milky Way might be detectable as they pass through 
detectors in laboratories on Earth. The very low cross section of WIMPs on ordinary 
material makes these interactions quite rare, but recent experiments have made progress. 
The direct detection experiments can measure and distinguish from background the tiny 
energy deposited by elastic scattering of a WIMP off a target nucleus. The current 
experimental results set limits on WIMP-nucleon cross sections, and we compare LNP- 
nucleon cross sections given in appendix [B] with these limits. Dark matter particles in 
the Milky Way halo have presumably a mean speed < v >~ 300 Km = 10~^c, 
therefore the process can be treated in the nonrelativistic limit. 

The nucleon coupling of a slow-moving Majorana fermion is characterized by two 
terms: spin-dependent (axial- vector) and spin-independent (scalar). We consider these 
two contributions separately. 

The spin-dependent cross section for LNP-nucleus elastic scattering is given by ()38p . 
For a proton target A'^J {J + 1) ~ 1 and the cross section is 



az {LN LN) 



3.5(^1/^20' X 10-^ pb 



(15) 



(for the definition of V see ([8])). The cross section (fTSl) for the three cases discussed 
above is always 2 — 3 orders of magnitude below current limits [T^. 

The spin-independent cross section is given by ()40p . It depends sensibly on the Higgs 
mass, so it is different between symmetric and asymmetric cases. For scattering from a 
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Figure 2: Spin-independent cross section: symmetric case (units 10~^p6) 
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(16) 



where ^ = V3; {aV2i + PVu). For the reason explained in section [5l we take the former 
reference value in the symmetric case and the latter in the asymmetric cases. The cross 
section (I16p for the symmetric case is plotted in figure [2] in units of 10"'' p5. It is about 
one order of magnitude above the experimental limits [13]. In the asymmetric cases, 
instead, spin-independent cross section is always 1 — 2 orders of magnitude below current 
limits, but within the sensitivity of experiments currently under study |14j . 

5 Higgs boson physics 

In this section we analyze the effects on Higgs boson physics induced by the new particles. 
In the first subsection we compute the contributions to the electroweak observables 
from their virtual exchanges and we will see that the upper limit on the Higgs mass 
is significantly affected. In the second subsection we analyze the new available decay 
channels to Higgs boson decays and compute the relevant branching ratios. 



5.1 Electro Weak Precision analysis 

The interaction lagrangian of the new particles with the gauge bosons is 
/\C\int = -Vu^ W+^-f^' Xi + h.c. + 
^24^/.^^7Vxi + h.c.+ 



2 
2 



i^ri'-^ iViiV2j + V2iVij) xiri' xi 



(17) 
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Figure 3: Region of the {S, T) plane allowed by EWPT at 68%CL and dependence of S 
and T on the Higgs mass. The thin black line marks ra^ = 400 GeV. Prom [15| . 

where sums under repeated indices are understood. The new particles contributions to 
T and S are respectively 

3 _ _ 

i=l 

1 ^ 

- 2 + V2^iVijf A K, -mj) (18) 

1 3 

5 = - ^ {VuV2j + V2iVi,f F {rrn, -mj) - F /i) (19) 

The functions F and A are defined in appendix O 

We have now all the ingredients to perform the analysis. We have verified that 
in the symmetric case the contribution to T is negligible, as required by the custodial 
symmetry mentioned in section [3l whereas in the asymmetric cases S is not significantly 
affected. The experimental contours in the (S, T) plane are shown in figure [3] and our 
results for the significant cases are shown in figure HI 

In the symmetric case AT is irrelevant and AS is inside the experimental ellipse for 
almost all the region that provides the entire dark matter abundance; if we raise the value 
of Yukawa couplings then AS" goes rapidly outside the ellipse, so we restrict ourselves 
just to this symmetric case and we do not consider higher values of a. Looking at figure 
[3] one can immediately see how an heavy Higgs can be allowed by ElectroWeak Precision 
Tests (EWPT): the only thing that we need is new physics producing a positive AT and 
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a not too large AS. To raise the Higgs mass up to 500 GeV the needed compensation 
is AT 0.2 [8j. The asymmetric case is perfectly suited to this purpose: it gives 
unimportant AS and a positive AT as desired. We first studied the case of Yukawa 
coupling A'^ equal to 1, then we raised its value until we reached AT ~ 0.2, and this 
corresponds to A*^ = 1.3 or equivalently a = 0.65 . All the values of A and A'^ that we 
consider are consistent with a Landau pole for the Yukawa coupling above the unification 
scale. When the Higgs boson mass is raised, however, the Higgs quartic coupling is very 
likely to have a Landau pole below the unification scale. 

5.2 Higgs boson decays 

Another important effect on the Higgs boson physics is the increase of its total width. 
In the parameter space region of interest to us the only new available decay channel is 
h XiXh since the decays into other new particles are kinematically forbidden. The 
partial width for such decay results in 

where the parameter ^ is defined in ([^^ . The Higgs total width predicted by the 
Standard Model T^*^ is known as a function of nih [E] . We consider four values for the 
Higgs mass, and the correspondent SM width are reported here 
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Figure 5: Branching ratio as a function of M for m/j = 120 GeV (left) and for rrih 
150 Ge\^ (right). Labels indicate the value of a. 
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The partial width r^^.^ for decay into two LNPs is given by (I20p . and thus we can compute 
the branching ratio 



BR{h^XX) 



XX 



-pSM _|_ p 
h 



(21) 



XX 



For nih = 120 GeV we compute the branching ratio in the symmetric case, whereas 
for higher masses we make the calculation for both the asymmetric cases. As seen from 
figure [Din both the asymmetric cases the only free parameter is M, since if we impose 
nih'^ = 0.105 the value of ^ is automatically fixed. This is not true for the symmetric 
case, where for each M there are up to three values of fx. The branching ratios are 
plotted as a function of M. In the asymmetric cases it is the only free parameter. In 
the symmetric case we consider the line of figure [T] corresponding to the lower value of 
the LNP mass. For the line corresponding to higher value of the LNP mass the decay is 
kinematically forbidden. For the line in the middle even small pole effects might modify 
the branching ratios considerably, because we are in a region where the phase space is 
nearly saturated. 

The branching ratios for rrih = 120 GeV and for nih = 150 GeV are plotted in figure 
[SJ In the simmetric case this channel dominates the total width. For higher values of 
the Higgs boson mass the branching ratios decrease, as a consequence that the Standard 
Model width increases faster than the partial width into two LNPs. For m/j = 200 GeV 
they are always below the 4%, while for higher values of the Higgs mass they are even 
smaller. 



6 Electric dipole moment 

We have taken the parameters (A, A*^, M, /x) to be real until now. We now explore the 
possibility of a CP violating phase. This phase could be present only in the symmetric 
case, since if one of the Yukawa couplings vanishes (as in the asymmetric case) all the 
parameters can be made real by a fields redefinition. In the general case we can redefine 
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Figure 6: Two loop contribution to the electric dipole moment of a fermion f 



fields so that (A, A'^, /x) are real, leaving a residual phase on the parameter M. The mass 
matrix found in ([6]) becomes 



Mn 






-\f2av 




(22) 



The phase Q induces an electron electric dipole moment (EDM) at two loops, the 
dominant diagram responsible for it is generated by charged and neutral particles and 
is shown in figure [6] ^ . The induced EDM moment is given by 



a^ruf 



■87r2< 



1=1 



m 



(23) 



w 



where 
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(24) 
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The matrix V diagonalizes the mass matrix and is such that V'^M^V = diag (mi, m2, ms) 
with real and positive diagonal elements. The sign on the right-hand side of equation 
()23p corresponds to the fermion / with weak isospin ±1/2 and / is its electroweak 
partner. 

The experimental limit on electron electric dipole moment at 95%CL level is [T7j 



141 < 1.7 X 10^ 



-27 



e cm 



(26) 
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Figure 7: Induced electron EDM for 6 = f > f i f > f v Green (lightest) regions are such 
that mi < 75 GeV, black shading indicates regions where the induced EDM is above 
experimental limit. 



We consider four different values of the phase, namely 6 = it/6, 7r/4, 7r/3, 7r/2, and we 
shade in the usual (/x, M) plane the regions where the induced EDM is above such limit. 
We identify in that plane also the region where the LNP mass is below 75 GeV, since this 
is the case of our interest. We also restrict ourselves to charged particle mass M below 
600 GeV as in the plot shown in figure [H The plots are shown in figure [71 The first 
result is that for M < 600 GeV we can have a LPN mass below 75 GeV only for small 
phases, otherwise the imaginary part not present before would require a cancellation in 
the mass matrix obtainable only for higher value of M. Regarding the induced EDM 
there are regions where it is above the experimental limit, but never inside regions such 
that mi < 75 GeV. On the contrary, for mi < 75 GeV, the induced EDM is always 
below the limit (I26p . It could however be accessible to next generation experiments [18] 

m- 



7 Conclusions 

In the last few decades it has been realized that the ordinary matter which we have 
been studying until now constitutes only about 5% of the total universe energy density. 
Evidence for nonluminous gravitating mass abounds on all the scales, from galactic 
to global ones of hundreds of Megaparsec. The measurements of the light element 
abundances and of the fluctuations in the Cosmic Microwave Background show that a 
significant part of the dark matter must be non baryonic. The Standard Model of particle 
physics does not contain such component. Another missing opportunity for the Standard 
Model is that gauge coupling unification does not occur at high energy. In this work we 
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have discussed a minimal extension of the Standard Model, which can explain all the 
observed dark matter abundance and improve the gauge coupling unification, focusing 
on the parameter space region for which the LNP mass is below mw ~ SOGeV. For 
such region the effects on the Higgs boson physics are worth of consideration. 

We have considered two limiting cases: almost equal Yukawa couplings (symmetric) 
and one of them vanishing (asymmetric). In both cases all the observed dark matter 
abundance could be explained by the LNP. We have computed also the full spectra of 
the model for all the cases, and they are consistent with negative searches from LEP2. 
The spin-independent direct detection cross section is above the current limits only for 
the symmetric case. In the asymmetric case it is well below these limits, as the spin- 
dependent cross sections for both cases. However they are all within the sensitivity of 
experiments currently under study. 

The new particles might have both direct and indirect effects on the Higgs boson 
physics. We have analyzed these effects and found that they are very different for each 
case. In the symmetric case the contribution to the electroweak observables is small, 
but the Higgs decays in LNP pairs dominate the total width. This might hide the 
Higgs boson at the Large Hadron Collider. On the contrary in the asymmetric case the 
contribution to the EWPT is important, and the indirect limit on the Higgs mass valid in 
the Standard Model can be raised. Finally we have considered a CP violating phase for 
the Dirac mass of the charge particle, giving rise to an electron electric dipole moment. 
We have verified that if we keep the LNP mass below 75 GeV the induced electric dipole 
moment is always below the current experimental limit, but perhaps accessible at the 
next generation experiments. 
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A Annihilations cross sections 

A.l s-channel Z exchange 

The LNP has a coupling with the Z boson given by 

Tp- {VuV2i) Z^Wn^j^'tPi (27) 

The cross section for the process Xl Xl ^* ^ f f for nonrelativistic LNPs and in the 
limit of massless final products is 

where the sum runs over all the Standard Model fermions except for the top quark. 
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A. 2 s-channel h exchange 

The coupling between the LNP and the Higgs boson is 

^Wi^ih ^ = V3i{aV2i+pVu) (29) 

The cross section for the process Xl Xl h* f f, kinematically identical to the 
previous one, is 

m^ml 2 ^2Q-j 

47r y2 (^4^2 _ ^2 ^ 2 rel 

where is the mass of the b quark (the process with bb in the final state is dominant, 
since ah oc rnj). 

A. 3 Thermally averaged cross sections 

To compute the thermally averaged cross section it is useful to observe that in both 
cases 

aj Vrel = bj vl^i (31) 

where j = Z, h, respectively for Z exchange and for h exchange. Performing the thermal 
average we obtain 

< ajVrel >= 6^ (32) 
so in our case the value of ctq is given by 

ao = 6{bz + bh) (33) 

B LNP-nucleus elastic cross sections 

B. l Spin-independent cross section 

The elastic scattering process is: xi M ^ Z* ^ xi ^ 

where TV is a generic nucleus. The vertex between quarks and the Higgs boson is given 
by the Standard Model Lagrangian and results in 

—Z^'Jl (34) 
cw 

The weak neutral current of the quarks J° is of the form X]^^7'^ (cy — c^7^) where 



the parameters cL e (?, are known as function of the Weinberg angle only. Wc have 



^1 

known as function of the Weir 
two different contributions to the amplitude: the quarks vector current and the quarks 
axial- vector current. We can neglect the first contribute for nonrelativistic LNP. The 
vector-axial contribution is described by the effective lagrangian 

Caxial = {V11V21) '^lYi^'^L ^1 ^^^''^ (35) 

Q 

where we define C = 2V2GFc'i. 
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We introduce the parameters 



A- [60 

The quantity J is the total angular momentum of the nucleus, < 5p > is the expectation 
value of the spin content of the proton group in the nucleus, and similarly for < Sn >■ 
The total cross section is 1201 



az {LN -^LN) = — {VuV2if GlmlK^J{J + 1) (38) 



vr 



where is the reduced mass of the system LNP-nucleus 

B.2 Spin-independent cross section 

The elastic scattering process is: xi ^ ^ h* ^ xi ^ 

The nucleonic matrix element can be parameterized by [21 



< iV| ^ mqqq\N >= f tun < N\N > / ^ 0.3 (39) 

and finally the spin-independent cross section results in 

„,(LN^LN) = ?^^ (40) 
where rrir is the reduced mass of the system LNP-nucleus. 



C ElectroWeak Precision Test 

As well known, new physics effects to the EWPT are conveniently represented by the 
parameters T and S, defined by 

^^ n33(o)-nw(o) s='^-^^i(,,{Q) (41) 

in terms of the vacuum polarization amplitudes 

mf/(Q) = r/^"% (g2) + q^^q- terms (42) 

with i, j = 3, 0,Ty for W^j", B'^, respectively. 

Expressions for the vacuum polarization amplitudes produced by fermions coupled 
to a generic gauge boson are known [22j. For a fermion loop with internal masses mi 
and 7712 and a vector coupling it is 



n(o) - ^ 



167r2 



(mi — m2f In —p^ — 2 ~ 2mim2 



mjmg 



2mim2 {m\ + m^ — mf — m^ m^ 
mf — mn mi 



(43) 
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n(o) 



247r2 



In 



mim2 (^3mf — 4mim2 + Sm^) 



+ 



In 



(44) 



For an axial coupling the results are obtained by letting mi — > —mi in the previous 
expressions. These results are valid for Dirac fermions, for Major ana fermions there is 
an extra factor of 2. A is a cutoff of the loop integral which disappears in the overall 
expressions (f^TI) . 

We define for convenience 

1 



A{mi,m2) 
F {mi, 1712) 



2a, 



,n(o) 



Ann (0) 



(45) 
(46) 
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